function [adjMatrix,nodeNames] = getDAG()
nodeNames{1} = 'HISTORY';
nodeNames{2} = 'CVP';
nodeNames{3} = 'PCWP';
nodeNames{4} = 'HYPOVOLEMIA';
nodeNames{5} = 'LVEDVOLUME';
nodeNames{6} = 'LVFAILURE';
nodeNames{7} = 'STROKEVOLUME';
nodeNames{8} = 'ERRLOWOUTPUT';
nodeNames{9} = 'HRBP';
nodeNames{10} = 'HREKG';
nodeNames{11} = 'ERRCAUTER';
nodeNames{12} = 'HRSAT';
nodeNames{13} = 'INSUFFANESTH';
nodeNames{14} = 'ANAPHYLAXIS';
nodeNames{15} = 'TPR';
nodeNames{16} = 'EXPCO2';
nodeNames{17} = 'KINKEDTUBE';
nodeNames{18} = 'MINVOL';
nodeNames{19} = 'FIO2';
nodeNames{20} = 'PVSAT';
nodeNames{21} = 'SAO2';
nodeNames{22} = 'PAP';
nodeNames{23} = 'PULMEMBOLUS';
nodeNames{24} = 'SHUNT';
nodeNames{25} = 'INTUBATION';
nodeNames{26} = 'PRESS';
nodeNames{27} = 'DISCONNECT';
nodeNames{28} = 'MINVOLSET';
nodeNames{29} = 'VENTMACH';
nodeNames{30} = 'VENTTUBE';
nodeNames{31} = 'VENTLUNG';
nodeNames{32} = 'VENTALV';
nodeNames{33} = 'ARTCO2';
nodeNames{34} = 'CATECHOL';
nodeNames{35} = 'HR';
nodeNames{36} = 'CO';
nodeNames{37} = 'BP';
adjMatrix = zeros(length(nodeNames));
adjMatrix(6,1) = 1;
adjMatrix(5,2) = 1;
adjMatrix(5,3) = 1;
adjMatrix(4,5) = 1;
adjMatrix(6,5) = 1;
adjMatrix(4,7) = 1;
adjMatrix(6,7) = 1;
adjMatrix(8,9) = 1;
adjMatrix(35,9) = 1;
adjMatrix(11,10) = 1;
adjMatrix(35,10) = 1;
adjMatrix(11,12) = 1;
adjMatrix(35,12) = 1;
adjMatrix(14,15) = 1;
adjMatrix(31,16) = 1;
adjMatrix(33,16) = 1;
adjMatrix(25,18) = 1;
adjMatrix(31,18) = 1;
adjMatrix(19,20) = 1;
adjMatrix(32,20) = 1;
adjMatrix(20,21) = 1;
adjMatrix(24,21) = 1;
adjMatrix(23,22) = 1;
adjMatrix(23,24) = 1;
adjMatrix(25,24) = 1;
adjMatrix(17,26) = 1;
adjMatrix(25,26) = 1;
adjMatrix(30,26) = 1;
adjMatrix(28,29) = 1;
adjMatrix(27,30) = 1;
adjMatrix(29,30) = 1;
adjMatrix(17,31) = 1;
adjMatrix(25,31) = 1;
adjMatrix(30,31) = 1;
adjMatrix(25,32) = 1;
adjMatrix(31,32) = 1;
adjMatrix(32,33) = 1;
adjMatrix(13,34) = 1;
adjMatrix(15,34) = 1;
adjMatrix(21,34) = 1;
adjMatrix(33,34) = 1;
adjMatrix(34,35) = 1;
adjMatrix(7,36) = 1;
adjMatrix(35,36) = 1;
adjMatrix(15,37) = 1;
adjMatrix(36,37) = 1;
P = topologicalPermutation(adjMatrix);
adjMatrix = adjMatrix(P,P);
nodeNames = nodeNames(P);
